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It is shown that the "hysteresis" in a polaron model of electron transport through the molecule 
found by M. Galperin et al. [Nano Lett. 5, 125 (2005)] is an artefact of their "mean-field" approx- 
imation. The reason is trivial: after illegitimate replacement n 2 = nno, where no = ^cjcoy ^ 1 
the average molecular level occupation Galperin et al obtained non-physical dependence of a non- 
degnerate molecular energy level on the non-integer mean occupation number no (i.e. the electron 
self-interaction) and the resulting non-linearity of current. The theory of correlated polaronic trans- 
port through molecular quantum dots (MQDs) that we proposed earlier [Phys. Rev. B67, 235312 
(2003)] proved that there is no hysteresis or switching in current-voltage characteristics of non- 
degenerate, d = 1, or double degenerate, d = 2, molecular bridges, contrary to the mean-field result. 
Switching could only appear in multiply degenerate MQDs with d > 2 due to electron correlations. 
Most of the molecular quantum dots are in the regime of weak coupling to the electrodes addressed 
in our formalism. 
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Although the correlated electron transport through 
mesoscopic systems with repulsive electron-electron in- 
teractions received considerable attention in the past, 
and continues to be the focus of intensive studies 0, 
much less has been known about a role of electron- 
phonon correlations in MQD. Recently we have pro- 
posed a negative— U Hubbard model of a d- fold degen- 
erate quantum dot and a polaron model of reso- 
nant tunneling through a molecule with d— degenerate 
levelQ. We found that the attractive electron correla- 
tions caused by any interaction within the molecule could 
provide a molecular switching effect where the current- 
voltage (I-V) characteristic has two branches with high 
and low current at the same voltage. This prediction has 
been confirmed by our theory of the correlated transport 
through degenerate MQDs with a full account of both 
the Coulomb repulsion and realistic electron-phonon (e- 
ph) interactions 0. We have shown that while the 
phonon side-bands significantly modify switching in com- 
parison with the negative-?/ Hubbard model (appearance 
of phonon ladder on the I-V curve), switching is robust. 
It shows up when the effective interaction of polarons is 
attractive and the state of the dot is multiply degenerate, 
d > 2, while there is no switching in a non-degenerate 
(d = 1) or a double degenerate (d — 2) MQD. 

Surprisingly, later on Galperin et al. 4] neglected these 
results, claiming that even a non-degenerate electronic 
level coupled to a single vibrational mode provides an I- 
V curve with the hysteresis, switching, and negative dif- 
ferential resistance. Here we show that these findings are 
artefacts of their mean-field approximation that neglects 
the Fermi-Dirac statistics of electrons. 



First, we will illustrate the failure of the mean-field 
approximation of Ref. Q on a simplest model of a single 
atomic level coupled with a single one-dimensional oscil- 
lator using the first quantization representation for its 
displacement x, 

„ . „ 1 ^ kx 2 
H = e a n + fxn-—— + — (1) 

Here M and k are the oscillator mass and the spring 
constant, / is the interaction force, and K = c = ks = 1- 
This Hamiltonian is readily diagonalized with the exact 
displacement transformation of the vibration coordinate 

x, 



x = y- hf/k, 



(2) 



to the transformed Hamiltonian without electron-phonon 
coupling, 



H 



e = 



en 
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2M dy 2 

£q — E p , 
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2 



(3) 
(4) 



where we used h 2 — n because of the Fermi-Dirac statis- 
tics. It describes a small polaron at the atomic level Eq 
shifted down by the polaron level shift E p = f 2 /2k, and 
entirely decoupled from ion vibrations. The ion vibrates 
near a new equilibrium position, shifted by f /k, with the 
"old" frequency (k/M) 1 / 2 . As a result of the local ion de- 
formation, the total energy of the whole system decreases 
by E p since a decrease of the electron energy by — 2E p 
overruns an increase of the deformation energy E p . It 
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becomes clear that the major error of the mean-field ap- 
proximation of Ref. 4] originates in illegitimate replace- 
ment of the square of the occupation number operator 
h = cjco by its "mean-field" expression n 2 = hno with 
the average population of a single molecular level, no, in 
disagreement with the exact identity, n 2 = h. This leads 
to a spurious self-interaction of a single polaron with it- 
self [e = £o — noEp instead of Eq.Q], and a resulting 
non-existent nonlinearity in the rate equation. 

The correct procedure should be as follows, see Ref.|3]. 
The appropriate molecular Hamiltonian includes the 
Coulomb repulsion, U c , and the electron- vibron inter- 
action as lal 



H = Ye h +-Y 



jjC - ~ 
U,,,,in n,j> 



+ IZ^^^A+- ff - c -)+II^(4 d 9 + 1 /2).(5) 

Here d q annihilates phonons, uj q is the phonon (vibron) 
frequency, and 7 are the e-ph coupling constant (q enu- 
merates the vibron modes). This Hamiltonian conserves 
the occupation numbers of molecular states h . 

One can apply the canonical unitary transforma- 



tion e s , with S 



E 



{l M d q ~ H - c ) integrating 



phonons out. The electron and phonon operators are 
transformed as 

5 M = c^X^, Xp = exp lw d i ~ H - c -j ( 6 ) 



and 



(7) 



respectively. This Lang-Firsov transformation shifts ions 
to new equilibrium positions with no effect on the phonon 
frequencies. The diagonalization is exact: 



where 



(8) 
(9) 



is the renormalized interaction of polarons comprising 
their interaction via molecular deformations (vibrons) 
and the original Coulomb repulsion, U£ ,. The molec- 
ular energy levels are shifted by the polaron level-shift 
due to a deformation created by the polaron, 



: -»-£bwp" 9 . 
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(10) 



Applying the same transformation to the retarded 
Green's function, one obtains the exact MQD spectral 



function Q for ad — fold degenerate MQD (i.e. the den- 
sity of molecular states, DOS) as 



d-l 



r=0 



2=0 



3 /W/2 [(! _ _ rU _ + n §^ - r u + l UQ )] 

+ (1 - S w )e- pulol/2 [nS(uj -rU - Ioj ) 

+ (1 - n)5{uo - rU + Iuj q )} , (11) 



where 



exp 



-|7| 2 coth 



(3uj 



(12) 



is the polaron narrowing factor |{|, £ = 
\j\ 2 / sinh(/3cj /2), Ii (0 the modified Bessel func- 
tion, [3 — 1/T, and 8^ the Kroneker symbol. To simplify 
our discussion, we assume that the Coulomb integrals 
do not depend on the orbital index, i.e. U^i = U, 
and consider a coupling to a single vibrational mode, 

UJq = UJ . 

The important feature of DOS, Eq.lJTTJ, is its non- 
linear dependence on the average electronic population 
n = (c^c^) which leads to the switching, hysteresis, and 
other nonlinear effects in I-V characteristics for d > 2 
3| • It appears due to correlations between different elec- 
tronic states via the correlation coefficients 



Z r (n) 



(d-l)! 
r\{d - 1-r) 



\d-l- 



(13) 



There is no nonlinearity if the dot is nondegenerate, 
d = 1, since Zo(n) = 1, contrary to Ref. Q. In this 
simple case the DOS, Eq. (|ll(l . is a linear function of 
the average population that can be found as a textbook 
example of an exactly solvable problems 0. As a result, 
the rate equation for n Q yields only a single solution 
see Eq. (43) in Ref. Q] and no switching in the whole 
voltage range. 

However, their "mean-field" approximation led the au- 
thors of Ref. 4] to the opposite conclusion. Indeed, 
Galperin et al. |4[ have replaced the occupation number 
operator h in the c-ph interaction by the average popula- 
tion no [Eq. (2) of Ref. 0] and found the average steady- 
state vibronic displacement (d + <$) proportional to no 
(this is an explicit neglect of all quantum fluctuations on 
the dot accounted for in the exact solution 0). Then, 
replacing the displacement operator d + $ in the bare 
Hamiltonian, Eq. (1), by its average, Galperin et al. ob- 
tained a new molecular level, io = £o — 2e reorg rto shifted 
linearly with the average population of the level. This 
is in stark disagreement with the conventional constant 
polaronic level shift, Eq. I|4I10[1 {e reorg is |7| 2 i^o in our no- 
tations). Their spectral function turned out to be highly 
nonlinear as a function of the population, e.g. for the 
weak-coupling with the leads p(lj) — S(LO—So—2e reorg no), 
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sec Eq. (17) inRef. Q. As a result, the authors of Ref. |4| 
have found multiple solutions for the steady-state popula- 
tion, Eq. (15) and Fig. 1, and switching, Fig. 4 of Ref. Q, 
which actually do not exist. Takinginto account the cou- 
pling with the leads (r in Ref. [3j,|^]) beyond the second 
order and the coupling between the molecular and bath 
phonons does not provide any non-linearity of the non- 
degenerate DOS either, because these couplings do not 
depend on the electron population. 

Note that the the mean-field solution by Galperin et 
al. applies at any ratio T/ujq, including the limit of 
interest to us, r -C ojq. where their transition between 
the states with uq = and 1 only sharpens, but none of 
the results change. Therefore, they do predict a current 
bistability in the system where it does not exist at d = 1. 
Results in Ref. 0] are plotted for T > u> , T w 0.1 - 0.3 
eV, which corresponds to molecular bridges with a resis- 
tance of about a few 100 kfL Such model 'molecules' are 
rather 'metallic' in their conductance and could hardly 
show any bistability because carriers do not have time to 
interact with vibrons on the molecule. This obvious con- 



clusion for molecules strongly coupled to the electrodes 
can be reached in many ways, see e.g. a very involved 
derivation in Ref. Q. The current hysteresis does not 
occur in their model, it remains a single-valued function 
of bias with superimposed noise. In any case, it cer- 
tainly has nothing to do with our result |3J that applies 
to molecular quantum dots (r -C u>o) with d > 2. 

As a matter of fact, most of the molecules are very 
resistive, so the actual molecular quantum dots are in 
the regime we study, see Ref.Q. For example the re- 
sistance of fully conjugated three-phenyl ring Tour-Reed 
molecules chemically bonded to metallic Au electrodes 
is larger than lGfi. Different from the non-degenerate 
dot, the rate equation for a multi-degenerate dot, d > 2, 
weakly coupled to the leads has multiple physical roots in 
a certain voltage range and a hysteretic behavior due to 
correlations between different electronic states of MQD 
13- We conclude that Galperin et al. have found 
a non-existent hysteresis in a model already solved well 
before their work in Ref. Q. 
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